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Abstract 

A challenging problem in aircraft engine health 
management (EHM) system development is to detect and 
isolate faults in system components (i.e., compressor, turbine), 
actuators, and sensors. Existing nonlinear EHM methods often 
deal with component faults, actuator faults, and sensor faults 
separately, which may potentially lead to incorrect diagnostic 
decisions and unnecessary maintenance. Therefore, it would 
be ideal to address sensor faults, actuator faults, and 
component faults under one unified framework. This paper 
presents a systematic and unified nonlinear adaptive 
framework for detecting and isolating sensor faults, actuator 
faults, and component faults for aircraft engines. The fault 
detection and isolation (FDI) architecture consists of a parallel 
bank of nonlinear adaptive estimators. Adaptive thresholds are 
appropriately designed such that, in the presence of a 
particular fault, all components of the residual generated by 
the adaptive estimator corresponding to the actual fault type 
remain below their thresholds. If the faults are sufficiently 
different, then at least one component of the residual 
generated by each remaining adaptive estimator should exceed 
its threshold. Therefore, based on the specific response of the 
residuals, sensor faults, actuator faults, and component faults 
can be isolated. The effectiveness of the approach was 
evaluated using the NASA C-MAPSS turbofan engine model, 
and simulation results are presented. 


Introduction 

National Transportation Safety Board accident data covering 
7,571 U.S.-registered aircraft from 1980 to 2001, categorized by 
accident cause, reveal that 52 percent of the hardware-induced 
accidents were related to aircraft system malfunctions and 
36 percent of these were caused by propulsion system 
component malfunctions (Ref. 1). Therefore, a real-time fault 
diagnosis scheme for aircraft engines might significantly 
improve flight safety by enabling accurate and early detection 
and isolation of incipient fault conditions. 

An important area of engine health management (EHM) is 
sensor validation. A sensor fault may lead to poor regulation 
or tracking performance, or even affect the stability of the 
control system. Moreover, a faulty sensor output may cause 
inaccurate diagnostics/prognostics, resulting in unnecessary 
replacement of system components or mission abortion. 
Therefore, it is important to correctly assess the health of on- 
board sensors. In addition to sensors, certain propulsion 
system components and actuators may fail as a result of aging 
or damage due to harsh operating conditions or combat. 
Existing nonlinear EHM methods often deal with component 
faults, actuator faults, and sensor faults separately. 
Specifically, when dealing with sensor validation, people 
usually assume there are no component or actuator faults; 
when dealing with component and actuator faults, it is often 
assumed that there are no sensor faults. In the former case, a 
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component or actuator fault may be misinterpreted as a sensor 
fault; in the latter case, a sensor fault may be misinterpreted as 
a component or actuator fault. Both cases may potentially lead 
to incorrect fault diagnostic decisions and unnecessary 
maintenance. Therefore, it would be ideal to address sensor 
faults, actuator faults, and component faults under one unified 
framework. Several researchers have investigated the 
development of such a unified fault diagnostic framework 
(Refs. 2 to 7). However, most of those results are based on 
linear engine models. 

In addition, the desire for future propulsion systems to 
perform over an extended range of operating conditions, 
characterized by dramatic variations in dynamic pressure and 
nonlinear thermal dynamics, requires enhanced performance 
and new functionalities. Consequently, the dynamics of 
aircraft engines are usually highly nonlinear and rapidly 
changing. Many existing fault detection and isolation (FDI) 
methods are based on the assumption that the system exhibits 
linear behavior in the neighborhood of a steady-state operating 
point, and therefore linearization-based methods are used (e.g., 
the Kalman filter-based methods developed by Merrill et al. 
(Ref. 8) and Kobayashi and Simon (Refs. 9 and 10). 
Approaches that involve linearizing the engine dynamics 
about several steady-state operating conditions and blending 
parameters and controllers for these operating points tend to 
be rather complicated. Moreover, when the effect of various 
faults is taken into account, the size and complexity of the 
scheduling and calibration tables are significantly increased, 
which makes design and real-time implementation very 
difficult. Therefore, future EHM designs will benefit 
significantly from new methods that are directly based on 
intrinsic nonlinearities of the engine dynamics. 

To address this issue, a unified nonlinear adaptive 
framework for detecting, isolating, and estimating sensor 
faults, actuator faults, and component faults for aircraft 
engines is developed in this paper. The presented approach is 
based on a bank of nonlinear fault diagnostic estimators that 
generate estimated measurements using advanced nonlinear 
adaptive estimation/learning techniques. The approach 
employs a nonlinear adaptive estimation architecture that is 
capable of directly dealing with nonlinear dynamic system 
models and nonlinear faults (Refs. 11). By extending the FDI 
logic to handle transient conditions as well as steady-state 
conditions, false alarms can be reduced, and early detection 
and isolation of sensor faults, actuator faults, and component 
faults can be achieved. 

With the presented adaptive nonlinear FDI techniques, 
unstructured modeling uncertainty is taken into account to 
improve performance over a wide range of operating regimes. 
Unstructured modeling uncertainty refers to the case where the 
modeling uncertainty function appears possibly in all state 
equations without being pre-multiplied by a known 
distribution matrix that satisfies certain conditions. This 
allows us to formally introduce adaptive thresholds for both 
the fault detection and isolation tasks. In general, adaptive 
thresholds have advantages over fixed thresholds because they 


enhance fault sensitivity and robustness with respect to 
modeling uncertainty (Ref. 15). The purpose of this paper is to 
highlight these features using a realistic gas turbine engine 
simulation and to demonstrate how sensor faults, actuator 
faults and component faults are dealt with under a unified 
framework. Preliminary simulation results are presented to 
show the effectiveness of the approach using a transient 
operating scenario. 

The remainder of the paper is organized as follows. First, 
the FDI approach is presented. Second, a neural network- 
based adaptive engine model is introduced to enhance fault 
sensitivity and diagnostic robustness with respect to modeling 
uncertainty resulting from normal engine degradation. Then, 
the application of the techniques to a high-fidelity engine 
simulation is presented. Simulation results for an example 
transient test case are provided. Finally, conclusions and 
possible future work are presented. 


Nomenclature 


C-MAPSS 

Commercial Modular Aero-Propulsion System 
Simulation 

D 

Modeling uncertainties in the output equation 

d 

Known bound for d 

EHM 

Engine Health Management 

FDE 

Fault detection estimator 

FDI 

Fault Detection and Isolation 

FIE 

Fault isolation estimator 

HPC 

High-pressure compressor 

HPT 

High-pressure turbine 

L 

The number of potential sensor faults 

LPC 

Low-pressure compressor 

LPT 

Low-pressure turbine 

N 

The number of potential component and actuator 
faults 

Nf 

Fan speed 

Nc 

Core speed 

NN 

Neural network 

P20 

Fan inlet pressure 

Ps30 

Combustor static pressure 

P50 

Core exhaust pressure 

SLS 

Sea-level static (standard day) 

T20 

Fan inlet temperature 

T24 

Low-pressure compressor discharge temperature 

T48 

HPT exhaust gas temperature 

T d 

Fault detection time 

TRA 

Throttle Resolver Angle 

U 

System inputs 

VBV 

Variable bleed valve 

vsv 

Variable stator vanes 

Wf 

Fuel flow rate 

X 

System state variables 

Y 

System outputs 

r/ 

Modeling uncertainties in the state equation 

?! 

Known bound for r| 
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Figure 1. — Structure of nonlinear FDI (Ref. 14). 


Nonlinear FDI Approach 

In this section, we give a brief description of the nonlinear 
FDI scheme (more details can be found in Refs. 11 to 14). As 
shown in Figure 1, the fault detection and isolation 
architecture is based on a bank of N + / + 1 nonlinear adaptive 
estimators, where N is the number of actuator and component 
faults under consideration, and / is the number of sensor faults. 
One of the nonlinear adaptive estimators is the fault detection 
estimator (FDE) used for detecting the occurrence of any 
faults, while the remaining N+l nonlinear adaptive estimators 
are fault isolation estimators (FIEs), which are activated for 
the purpose of fault isolation only after a fault is detected. 

Consider a generic nonlinear jet engine system described 
by: 

x — Ax + y(x, u) + tj(x, u , t) + p x (t - T x )^(x, u) 
y = Cx + d(x,u,t)+ /3 y (t -Ty)F0 ^ 

where x, u , and y are the state vector, input vector, and output 
vector of the engine system, respectively, ( A,Q is an 
observable matrix pair, and t represents time. The model 

X„ = Ax n +Mx n ,u) 
y„ = Cx n 

represents the known nominal nonlinear dynamics of a typical 
new engine derived based on engine dynamics, while the 
healthy normal engine system is modeled by 

X-h = Ax h + A(x h ,u)+Jl(x h ,U,t) 

yh = Cxfj + d(xf,,u,t) ^ 


The differences between the nominal model of Equation (2) 
and the healthy system model of Equation (3) are the terms r\ 
and d , which represent the modeling uncertainties in the state 
equation and output equation, respectively. For instance, rj 
captures the changes in the engine dynamics as a result of 
normal degradation, and d includes disturbances in the sensor 
signal. 

The change in the system dynamics as a result of a sensor 
fault is characterized by J3 y {t - T y )F6 in Equation (1). 
Specifically, the vector FO represents a sensor bias fault, and 
the function J3 y (t - T y ) characterizes the time profile of the 
sensor fault, where T y is the unknown fault occurrence time. 
The class of sensor bias faults under consideration is 
represented by FO, where F is the fault distribution vector, and 
6 is the unknown magnitude of the sensor bias. For instance, 
in the presence of a single sensor fault, the sensor fault 
distribution vector F has only one non-zero entry, which 
represents the corresponding corrupted output measurement. It 
is assumed that there are / sensors being monitored. 
Depending on the location of the fault, the distribution vector 
F belongs to a class of possible vectors {F l ,F 2 ,..., F 1 }, where, 
for any j = 1 ,...,/, only the yth component of the vector F j takes 
the value of 1, while all other components are zero. 
Accordingly, the scalar 0 is the unknown bias in the 
corresponding sensor. 

The change in the system dynamics as a result of a 
component fault or actuator fault is characterized by 
Pxif- T x )(/j{x,u) in Equation (1). Specifically, J3 x (t- T x ) denotes 
the time profile of the component or actuator fault that occurs 
at some unknown time T x , and <Kx,u) represents the nonlinear 
component or actuator fault function, which is modeled as a 
nonlinear function of x and u. It is assumed that there are N 
types of possible actuator and component faults in the fault 
class. Each fault function f\p= is described by 
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i/>p(x,u)= [6\P ) T gf (x,u),...,(0 m pf gP(x,u) 


( 4 ) 


where Of , i = 1 ,...,m, is an estimated parameter vector 
assumed to belong to a known compact and convex set, gf is 
a known smooth vector field, and m is the number of system 
state variables. As discussed in (Refs. 11 to 14), the fault 
model described by Equation (4) characterizes a general class 
of nonlinear faults where the vector field gf represents the 
functional structure of the pth fault affecting the zth state 
equation, while the unknown parameter vector Of 
characterizes the corresponding unknown fault magnitude. 

In this research, the following assumptions are made: 

1) Only one fault occurs at a time; 

2) Only the case of abrupt (sudden) faults is considered. 
Therefore, J3 x (t - T x ) and fi y (t - T y ) each take the form of a 
step function; 

3) The effect of normal engine degradation is considered 
as modeling uncertainty, not a fault; 

4) Modeling uncertainties, represented by rj and d in 
Equation (1), are unstructured and unknown nonlinear 
functions of x, u , and t , but bounded, i.e., 
rj(x,u,t)<rf(y,u,t), and d(x,u,t)< d(y,u,t ) . The 
known bounds rj and d on the unstructured modeling 
uncertainties are needed in order to distinguish between 
the effects of faults and modeling uncertainties 
(Refs. 11 to 14). 


functional structure of a particular fault type. The adaptation 
in the FIEs arises due to the unknown fault magnitude (i.e., 0 
in Equation (1) and Of in Equation (4)). The adaptive law for 
adjusting the fault magnitude estimate is derived using a 
Lyapunov synthesis approach (see, for instance, (Refs. 16 to 
18). Let us denote V s (t) as the estimated output vector 
generated by the sth FIE, where s e {1, — , N + 1} . The 
adaptive threshold for each FIE is designed in such a way that 
in the presence of the sth fault, the residual generated by each 
FIE may exceed its corresponding threshold, except for that of 
the sth FIE. 

The fault isolation scheme is developed following the 
general method presented in (Refs. 1 1 to 14). Each of the N + l 
isolable faults will have an FIE for which a set of / adaptive 
thresholds can be designed. Thus, if fault s occurs and is 
detected at time instant T d , then all / components of the 
residual vector from the corresponding FIE (the FIE with the 
fault s hypothesis) will remain below their adaptive thresholds 
for all t > T d . Likewise, if at least one component of the 
residual vector of FIE s exceeds its adaptive threshold for 
some finite time t > T d , then the possibility of the occurrence 
of fault s can be excluded. In other words, the FDI logic 
indicates that “the sth fault is isolated” if the following two 
conditions are met: 1) All the diagnostic residuals generated 
by FIE s remain below their corresponding adaptive 
thresholds; and 2) At least one component of the diagnostic 
residual generated by each remaining FIE exceeds its 
corresponding threshold. 


Fault Detection Method 

The design of the fault detection estimator is based on 
nonlinear adaptive estimation techniques (Refs. 11 to 14). The 
FDE continuously monitors the behavior of the engine system 
based on the nominal dynamical model given by Equation (2). 
The FDE is designed with the hypothesis that no fault is 
present in the system, so if at least one of the measured 
variables of the real system significantly deviates from the 
corresponding estimate, then fault occurrence can be 
concluded. Therefore, the fault detection time T d is the time 
when a residual first exceeds its threshold. 

If a fault is detected, it can be a sensor fault, a component 
fault, an actuator fault, or even some other anomalous 
condition. Isolation of faults (i.e., determining if a fault in the 
fault class under consideration has occurred) is accomplished 
by the fault isolation procedure described below. 

Fault Isolation Method 

As discussed above, the faults under consideration include / 
types of sensor faults and N types of actuator and component 
faults. Therefore, there are N + / faults in the augmented fault 
class. The bank of FIEs is activated only after a fault is 
detected. Each adaptive FIE is designed based on the 


Adaptive Thresholds for Fault Detection and 
Isolation 

It is well-known that modeling uncertainty has a corruptive 
effect on any model-based diagnostic method. As described 
earlier, the terms rj and d in Equation (3) represent the 
modeling uncertainties in the state equation and output 
equation, respectively. In aircraft engine models, uncertainties 
include the effect of normal engine health degradation, engine- 
to-engine variation, parameter variations due to dynamic 
operating environments, and sensor noise. Due to the presence 
of such modeling uncertainties, residuals are never zero even 
in the absence of faults. Moreover, the residual may change 
with the varying control inputs and dynamic operating 
conditions of the engine. Therefore, using a small fixed 
threshold may result in significant false alarms, while using a 
large fixed threshold may increase the number of missed 
detections/isolations, as illustrated in Figure 2. 

In this nonlinear FDI methodology, adaptive thresholds in 
the form of dynamic nonlinear functions of the system inputs, 
outputs, and estimated parameters are rigorously established 
based on the bounding functions rj and d on the modeling 
uncertainties and adaptive estimation techniques (Refs. 1 1 to 
14). The FDI thresholds automatically adapt to the changing 
engine operating conditions and engine dynamics to enhance 
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Figure 2. — Concept of adaptive thresholds in FDI. 


control inputs 



sensor outputs 



Figure 3. — Structure of the adaptive reference engine model. 


the robustness and fault sensitivity of the FDI scheme. The 
adaptive thresholds characterize the effect of several sources 
of uncertainty on the FDI residual, such as modeling 
uncertainty 77 , measurement uncertainty d , and the estimation 

error of fault parameters (i.e., ^in Eq. (1) and Of in Eq. (4)). 
Intuitively, the smaller the uncertainty is (resulting in a 
smaller threshold), the easier the task of isolating the 
associated faults. On the other hand, the capability to isolate a 
fault depends not only on the significance of uncertainty, but 
also on the degree that the faults are mutually “different” from 
each other, characterized by the so-called fault mismatch 
function (Refs. 11 to 14). 

Enhancement of Robustness and Fault 
Sensitivity 

In this research work we also use an adaptive nonlinear 
reference engine model to improve the accuracy and robustness 
of the diagnostic method. 


Similar to faults, normal engine health degradation causes the 
engine outputs to deviate from their reference values, which 
may be misinterpreted as a fault and consequently generate false 
alarms (Ref. 10). To compensate for deviations from the 
nominal engine model due to normal aging and engine-to- 
engine variation, this FDI method uses an adaptive on-board 
nonlinear engine reference model as shown in Figure 3. 
Specifically, the reference engine model consists of a nominal 
engine model, representing the dynamics of a typical new 
engine, and a neural network- (NN-) based adaptive component 
for capturing various modeling errors. 

One important issue that limits the practicality of adaptive 
reference engine models is the complexity involved in the 
update process. In the present method, the update process is 
rather simple. Specifically, the nominal engine model is fixed 
once designed, and weights of the neural network component 
are initialized to zero. Nonlinear adaptive approximation 
techniques (Ref. 16) are used to train the neural network on-line 
during flight to capture changes in engine dynamics as a result 
of normal aging (represented by p(x,uf) and d(x,u,t) in Eq. (1)). 
If no fault is detected after a number of flights, then the 
reference engine model can be simply updated off-line by 
replacing the old neural network weights in the adaptive engine 
model with the new ones recently learned. Hence, the reference 
engine model is always a good representation of the current 
engine dynamics. To summarize, the learning algorithm of the 
neural network runs in real-time during flights, but the reference 
engine model is updated off-line. 


Application to a Commercial Engine 
Platform 

In this section, the generic fault diagnostic techniques 
described above are applied to a practical engine fault 
diagnosis problem. The NASA C-MAPSS engine model 
(Ref. 19) was utilized as an example to evaluate the 
performance of the presented methodology. This evaluation 
was accomplished by incorporating the FDI code into the 
C-MAPSS platform to conduct the simulation evaluations. 

The C-MAPSS engine is a realistic representation of a 
90,000-lb thrust class turbofan engine with high-bypass ratio 
and a two-spool configuration. The simulation operates with a 
representative closed-loop controller that regulates fan speed 
(Nf) by manipulating three actuators (fuel metering valve, 
VSV and VBV). Figure 4 shows the various components of 
the engine. In this engine simulation, the five rotating 
components consist of the fan, low pressure compressor(LPC), 
and low pressure turbine (LPT) on the fan shaft, and the high 
pressure compressor (HPC) and high pressure turbine (HPT) 
on the core shaft. The station designations are shown at the 
bottom of the figure. As shown in Table 1, C-MAPSS 
provides a group of health parameters that can be adjusted to 
simulate damage, deterioration, or fault conditions. The 
dynamics of several actuators and sensors have also been 
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Figure 4. — Diagram of the C-MAPSS engine components and station designations (Ref. 20). 


included in the simulation (see DeCastro et al. (Ref. 20)). 
Table 2 shows the complete set of actuators and sensors 
included in the simulation and their respective bandwidths. 


TABLE 1.— HEALTH PARAMETERS 



Parameter 

Description 

1 

mod, Fan 

Fan efficiency modifier 

2 

^mod,Fan 

Fan corrected flow modifier 

3 

R^mod,Fan 

Fan pressure ratio modifier 

4 

^7 mod, LPC 

LPC efficiency modifier 

5 

^mod,LPC 

LPC corrected flow modifier 

6 

R^mod,LPC 

LPC pressure ratio modifier 

7 

V mod,HPC 

HPC efficiency modifier 

8 

^mod,HPC 

HPC corrected flow modifier 

9 

R^mod,HPC 

HPC pressure ratio modifier 

10 

V mod, HPT 

HPT efficiency modifier 

11 

^mod,HPT 

HPT corrected flow modifier 

12 

?7mod,LPT 

LPT efficiency modifier 

13 

^mod,LPT 

LPT corrected flow modifier 


TABLE 2.— ACTUATOR AND SENSOR 
DESCRIPTIONS AND BANDWIDTHS 


Variable 

Description 

Bandwidth 

Wf 

Fuel metering valve 

25 rad/s 

VBV 

Variable bleed valve 

23 rad/s 

VSV 

Variable stator vanes 

23 rad/s 

Nf 

Fan speed 

50 rad/s 

Nc 

Core speed 

50 rad/s 

T20 

Fan inlet temperature 

9 rad/s 

P20 

Fan inlet pressure 

25 rad/s 

T24 

LPC outlet temperature 

9 rad/s 

Ps30 

Combustor static pressure 

25 rad/s 

T48 

HPT exhaust gas temp. 

9 rad/s 

P50 

Core exhaust pressure 

25 rad/s 
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The parameters in Table 3 are used in our implementation 
to simulate sensor noise. 


TABLE 3.— PARAMETERS FOR THE SENSOR 
NOISE AND DYNAMICS MODULES 


Sensor 

Standard deviation 

SLS* max power 

Nf 

- 

2388 rpm 

Nc 

- 

9048 rpm 

T20 

1.0 °R 

519 °R 

T24 

2.2 °R 

642 °R 

T48 

2.2 °R 

2070 °R 

P20 

0.032 psi 

14.6 psi 

Ps30 

2.8 psi 

523 psi 

P50 

0.28 psi 

19 psi 


*SLS - Sea Level Static 


Computation of the Adaptive Thresholds 

As the engine experiences degradation with usage, the 
modeling discrepancies can cause the diagnostic performance 
to degrade; therefore adaptive FDI thresholds are employed to 
enhance diagnostic performance. Computation of the adaptive 

thresholds requires functions rj and d that bound the 

uncertainties in the state equation and output equation, 
respectively. Since analytical forms of these equations do not 
exist in numerical engine simulations such as C-MAPSS, a 
more direct modeling approach must be employed. In the 
modeling scheme shown in Figure 5, the healthy engine model 
is an adaptive model that includes a NN component 
representing normal engine degradation. The degraded 
C-MAPSS engine model is fed by an assumed “worst-case” 
health parameter combination (Ah) that can occur during the 



NN training cycle. As such, the functions 7j and d are taken 
as representations of the worst-case mismatch (uncertainty) 
between the model and actual engine. These dynamic 
functions are characterized by the degradation incurred during 
a small number of flights and are utilized to calculate the 
adaptive thresholds at each time step. 

Engine FDI Algorithm 

A proof-of-concept FDI algorithm, designed to detect and 
isolate six faults, was implemented and integrated into the C- 
MAPSS engine model. The six fault types are listed in 
Table 4. As shown in Figure 6, the prototype software consists 
of one FDE and six FIEs, one for each fault type. The sensor, 
actuator, and component faults considered in this work are a 
representative subset of faults that are considered relatively 
challenging to detect/isolate. For instance, the Nf sensor fault 
mode was selected because the primary engine fuel flow (Wf) 
controller uses the Nf variable for feedback. Therefore, Nf 
sensor faults may be masked by the control system as it 
regulates that variable. The sensed Nc and T24 values are used 
for scheduling the variable stator vanes (VSV). Therefore, 
studying the Nc and T24 sensors is of interest because shifts in 
those sensors can be possibly misinterpreted as a fault in VSV 
(and vice versa). Lastly, the HPC component fault was 
selected because it can possibly be misinterpreted as a fault in 
the VBV or VSV actuators due to the physical coupling. In 
this paper, HPC faults are modeled simply as commensurate 
step decreases in the HPC efficiency and flow capacity. 

The engine FDI logic is as follows. The engine FDI 
declares that “fault s is isolated” if the following two 
conditions are both met: 

• All the diagnostic residuals generated by FIE s remain 
below their corresponding adaptive thresholds. 

• At least one component of the diagnostic residual 
generated by each remaining FIE exceeds its 
corresponding adaptive threshold. 

Note that the fault can be a sensor fault (Nf, Nc, and T24), 
an actuator fault (Wf and VSV), or a component fault (HPC). 
This FDI logic is illustrated schematically in Figure 7. Note 
that in this figure three residual signals (Nf, Nc and T24) are 
generated by each FIE as an example. The number of residual 
signals may be increased when more sensor faults need to be 
considered. Particularly, “0” indicates the residual always 
remains below its threshold, while “1” indicates the residual 
exceeds its threshold at some finite time. When there is only 
one “0” in the isolation indicator, the fault is successfully 



Figure 5. — Approach for determining parameters used to 
compute adaptive thresholds. 



Figure 6. — A prototype FDI implemented for proof-of- 
concept demonstration. 


TABLE 4. — FAULTS CONSIDERED 
IN THE SIMULATION STUDY 


Fault type 

Failure mode 

Nf sensor 
Nc sensor 
T24 sensor 
Wf actuator 
VSV actuator 
HPC 

Bias in sensor reading 

Bias in sensor reading 

Bias in sensor reading 

Bias in actuator output 

Bias in actuator output 

Abrupt fault resulting in changes of the 

health parameter values in HPC 
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isolated. If there is no “0”, the fault detected is an unexpected 
fault. When there are multiple “0”s, the fault is not isolable 
among the faults whose corresponding FIEs give “0”s. Note 
that the fault isolability depends on the interplay between 
modeling uncertainty and the difference between faults. The 
fault isolability conditions have been rigorously investigated 
by using the so-called fault mismatch function, which 
characterizes the differences between faults (Refs. 11 to 14). 


FIE1 (Nf) 
FIE2 (Nc) 
FIE3 (T24) 
FIE4 (Wf) 
FIE5 (VSV) 
FIE6 (HPC) 



Figure 7. — Implementation of fault isolation logic. 
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Figure 8. — Time history of the transient scenario 
considered. 


Simulation Results — Transient Example 

To evaluate the proof-of-concept nonlinear FDI scheme, a 
simulation study for all six fault modes was successfully 
conducted and demonstrated using the C-MAPSS engine 
model. In this initial work, only one transient example was 
explored. To test the approach in a time- varying operating 
scenario, a transient scenario similar to that in Kobayashi and 
Simon (Ref. 10) was used. The throttle resolver angle (TRA) 
was ramped from 60° to 69° in 1 sec. After steady-state 
operation, the TRA was ramped back to 60° in 1 sec, as shown 
in Figure 8. The altitude and Mach number were held constant 
at 25000 ft and 0.62, respectively. In this example, the 
difference between the engine model and the real engine (for 
instance, due to normal degradation) is a represented by a 
-0.2 percent shift in each of LPC efficiency, LPC flow 
capacity, HPC efficiency, and HPC flow capacity. 

Simulation results were generated for all six faults with 
different fault sizes. The FDI software successfully detected 
and isolated all six faults as long as the simulated fault 
magnitude was larger than a minimum detectable/isolable size 
(to be discussed next), and no misdiagnosis occurred. In 
addition to fault diagnosis, the adaptive FIEs also provide an 
estimate of the fault size (0 in Eq. (1) and Of in Eq. (4)), 
which can potentially be utilized for automated fault 
accommodation. 

For the sake of space, only the simulation results of the Nf 
(fan speed) sensor bias fault are presented here as an 
illustrative example. In this simulation study, the neural 
network-based adaptive element is not activated. Figures 9, 10 
and 11 correspond to the case of a 20 rpm bias (nearly 
1 percent of the simulated operating condition) in the Nf 
sensor signal, which occurred at t = 35 sec. As can be seen in 
Figure 9, the fault was detected at t = 35.1 sec. The estimated 
Nf fault magnitude is 22 rpm (shown in Fig. 10), which is very 
close to the actual fault size (20 rpm). The small discrepancy 
in the sensor bias estimation is due to the modeling error 
simulated by introducing nominal engine degradation. The 
steady-state bias estimate is even closer to the actual fault size 
when no modeling error is present. Figure 11 shows the time 
history plots of the three residuals generated by each of the six 
FIE’s, along with their corresponding residual thresholds. 
Since the FIEs are activated only if a fault is detected by the 
FDE, the time scale on all six subplots begins at 1=35.1 sec. As 
shown in Figure 11, the single FIE for which all the residuals 
always remain below their corresponding thresholds was 
designed on the hypothesis of an Nf bias fault, while at least 
one component of the diagnostic residuals generated by each 
remaining FIE exceeds its threshold almost immediately 
following the positive fault detection decision. By applying 
the fault isolation logic shown in Figure 7, an Nf bias is 
isolated by the FDI algorithm. 
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FDE residuals 




Figure 9. — Residuals generated by the FDE and adaptive 
thresholds (Nf bias = 20, fault occurrence time = 35 sec). 


In addition, a study of the minimum detectable and isolable 
fault size for all six fault modes was conducted. The results 
are summarized in Table 5, using the transient scenario in 
Figure 8 as before. From this study, the following observations 
can be seen. First, the results indicate that the nonlinear FDI 
can correctly detect and isolate faults in each of the chosen 
sensor, actuator, and engine components. Second, in each 
case, the minimum detectable and isolable faults are within 
5 percent of the steady-state value and many are significantly 
less, which is a promising result for detecting small-scale 
though process noise and a transient operating scenario were 
considered in this evaluation. Initial development of a NN- 
based adaptive reference engine model was also conducted for 



Figure 10. — Nf sensor fault bias estimation. 


faults. Third, in no case was a false alarm generated even 
enhancing fault sensitivity and robustness to modeling 
uncertainty. As an example, a preliminary result is shown in 
Figure 12. In this example, the simulated engine degradation 
includes 0.3 percent losses in fan efficiency and flow capacity, 
0.3 percent losses in LPC efficiency and flow capacity, and 
0.3 percent losses in HPC efficiency and flow capacity. No 
faults were present in the system. The NN was trained during 
the first 15 sec. The training was then stopped, and then the 
NN with the trained weights was utilized for the rest of the 
simulation. In Figure 12, it is obvious that the Nc residual was 
reduced significantly after the NN was trained and utilized in 
the engine model. 

The NN can be trained either off line or online in parallel 
with the FDI system but is not used in the FDI system while 
the NN is being trained. The NN in the FDI system is updated 
as necessary after a number of flights. The adaptive thresholds 
are affected when the adaptive engine model is updated. 


TABLE 5.— MINIMUM DETECTABLE AND ISOLABLE FAULT SIZE 


Fault type 

Minimum detectable fault size 
(% of SLS) 

Minimum isolable fault size 
(% of SLS) 

Estimated fault size 

Nf sensor bias 

5 (0.21%) 

10 (0.42%) 

11.2 

Nc sensor bias 

16 (0.18%) 

25 (0.28%) 

23.7 

T24 sensor bias 

5 (0.78%) 

11 (1.6%) 

11.1 

Wf actuator bias 

-120 (0.49%) 

-125 (0.52%) 

-110 

VSV actuator bias 

-0.20 (3.39%) 

-0.25 (4.24%) 

-0.21 

HPC efficiency and flow capacity loss 

-0.009 (0.9%) 

-0.015 (1.5%) 

-0.016 
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T24 Residual Nc Residual Nf Residual T24 Residual Nc Residual Nf Residual T24 Residual Nc Residual Nf Residual 
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(c) FIE no. 3; the Nf residual exceeds threshold 
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(b) FIE no. 2; the Nf residual exceeds threshold 
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(d) FIE no. 4; the Nf residual exceeds threshold 
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(e) FIE no. 5; the Nf residual exceeds threshold 
Figure 1 1 . — Sensor residuals generated by the FIEs 
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(f) FIE no. 6; the Nf residual exceeds threshold 
and their adaptive thresholds for Nf sensor fault scenario. 
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Figure 12. — NN enhancement simulation results. 

Conclusions and Future Work 

In this paper a new methodology for model-based fault 
detection and isolation (FDI) that is based on a unified 
nonlinear adaptive framework is described. The FDI 
architecture is capable of detecting and isolating sensor faults, 
actuator faults, and component faults in aircraft engines using 
a parallel bank of nonlinear adaptive estimators. The 
framework makes use of adaptive thresholds that intrinsically 
capture the nonlinear behavior of the engine, thereby 
addressing the limitations of fixed or scheduled thresholds. 
The approach was evaluated on the C-MAPSS engine model, 
revealing the promise of the approach for engine FDI. 

Although the results are promising, the evaluation of the 
algorithm presented in this paper is preliminary. In order to 
gain confidence in the technique, it is necessary to perform the 
following tasks. First, to compensate for modeling errors such 
as normal aging and engine-to-engine variation, it will be 
necessary to develop an adaptive real-time nonlinear engine 
model using neural network (NN) techniques, as discussed in 
this paper. The preliminary development of NN-based 
adaptive reference engine models will need to be further 
explored. Second, it is necessary to redevelop the FDI scheme 
by considering the full suite of sensors, actuators, and engine 
component faults. Once the full FDI system is developed, it is 
necessary to statistically characterize the performance of the 
algorithm at various flight conditions and levels of 
degradation and evaluate the false alarms and missed 
detections. Lastly, deployment of the code on a physical test 
bed is necessary for maturation of the technology. 
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